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Abstract 

In this paper we investigate the general Ericksen-Leslie (E-L) system modeling nematic 
liquid crystal flows. First, by overcoming the difficulties from lack of maximum principle for 
the director equation and high order nonlincarity for the stress tensor, we prove existence of 
global-in-time weak solutions under physically meaningful boundary conditions and suitable 
assumptions on the Leslie coefficients, which ensures that the total energy of the E-L system 
is dissipated. Moreover, for the E-L system with periodic boundary conditions, we prove the 
local well-poscdncss of classical solutions under the so called Parodi's relation and establish 
a blow-up criterion in terms of the temporal integral of both the maximum norm of the curl 
of the velocity field and the maximum norm of the gradient of liquid crystal director field. 
Key words: Nematic liquid crystal flow, Ericksen-Leslie system, existence of weak solutions, 
blow up criterion. 
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1 Introduction 

The hydrodynamic theory of liquid crystals due to Ericksen and Leslie was developed around 
1960's [5]IB]118| . The Ericksen-Leslie system for liquid crystal flows was derived from the macro- 
scopic point of view and was very successful in understanding the coupling between the velocity 
field and the molecular director field, especially in liquid crystals of nematic type. Assume 
that the material occupies a bounded spatial domain C M 3 with smooth boundary T. We 
denote by u = (ux,U2,us) T the velocity field of the flow and by d = (di,d2,d^) T the director 
field, which stands for the averaged macroscopic/continuum molecular orientation in M 3 . The 
resulting PDE system that we are going to consider in Qt '■= (0, T) x £1 can be written as [22|l41j : 

u t + (u • V)u + VP = -V-(VdOVd) + V-er, (1.1) 
V-u = 0, (1.2) 

d t + (u • V)d - cod + ^Ad = ~l(Ad- V d W(d)). (1.3) 
Ai Ai 

Equations (jl.ip and (|1 .3j) represent the conservation of linear momentum and angular momen- 
tum, respectively. Here, we consider the flow of an incompressible material that satisfies the 
incompressiblility condition (j 1 . 2 j) . P is a scalar function representing the hydrodynamic pres- 
sure (including the hydrostatic part and the induced elastic part from the director field). The 
term — Vd©Vd is the elastic (Ericksen) stress tensor, where Vd©Vd denotes the 3x3 matrix 
whose (i,j)-th entry is given by Vjd • Vjd, for 1 < i,j < 3. The symbol V r indicates the 
transpose of the gradient and 

A = ^(Vu + V T u), w = ^(Vu-V T u) 

represent the rate of strain tensor and the skew-symmetric part of the strain rate, respectively. 
For the sake of simplicity, we denote by 

d = d t + (u • V)d, N = d - u d = -^Ad - ^ (Ad - V d W(d)) (1.4) 

Ai Ai 

the material derivative of d (transport of center of mass) and the rigid rotation part of the 
changing rate of the director by fluid vorticity. The kinematic transport X\M + A2^4d = Ad — 
V c jW r (d) for the molecule director represents the effect of the macroscopic flow field on the 
microscopic structure. The material coefficients Ai and A2 reflect the molecular shape (Jeffreys 
orbit) and how slipper the particles are in the fluid (cf. |14^33j). The viscous (Leslie) stress 
tensor cr has the following form [4^119] : 

a = m(d T Ad)d © d + ^ © d + /i 3 d <g) N + n^A + ^{Ad) (g> d + fi 6 d © (Ad), (1.5) 

where (g> stands for the usual Kronecker product, i.e., (a ® h)ij := aibj , for 1 < i,j < 3. 
The six independent coefficients fix, fj,Q , which may depend on material and temperature, are 
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called Leslie coefficients. These coefficients are related to certain local correlations in the fluid 
(cf. e.g., [3] for more details). The function W penalizes the deviation of the length |d| from 
the value 1, which is due to liquid crystal molecules being of similar size (cf. [2TJE21E7] ) • The 
notation Vd represents the gradient with respect to the variable d . A typical example of W is 
the Ginzburg-Landau approximation that has a double well structure |21[I22]: 

^(d) = i L(|d| 2 -l) 2 . (1.6) 

The system (|l.ip — (|1.3p should be supplemented with suitable initial and boundary conditions. 
Here, we assume that it is subject to the following initial conditions 

u\ t =o = u with V • u = 0, d| t= o = d , in O, (1.7) 

and the homogeneous Dirichlet boundary condition for the velocity field 

u = 0, on(0,T)xr, (1.8) 

together with the nonhomogeneous Dirichlet boundary condition for d (constant in time for 
simplicity) 

d| r = d |r, on(0,T)xr. (1.9) 

The reformulated E-L system (|l.ip - (|1.3p with Ginzburg-Landau approximation W was in- 
troduced in [22] motivated by the work on the harmonic heat flow (see, e.g., |21[ I31] for detailed 
analysis on a simplest liquid crystal system, in which the Leslie stress a was neglected except the 
fluid viscosity term involving /Z4). It can also be derived via an energetic variational approach 
(cf. [2]), based on the basic energy law of the system under proper assumptions on the Leslie 
coefficients (e.g., the Parodi's relation, see (|2.5p below). Due to the complicated mathematical 
structure of the E-L system, most of the previous works were restricted to certain simplified ver- 
sions, see for instance, [2^ fTTj[2"Tl -[2"3 , 25 , 26 , 28 , 31 , 33 , 40j and the references therein. In particular, 
existence of weak solutions to the E-L system (|l.ip - (|1.3|) (with Ginzburg-Landau approxima- 
tion) subject to boundary conditions of Dirichlet type (|1.8|) — (|1.9|) was obtained in [22] under 
the specific assumption that A2 = 0. Physically this assumption indicates that the stretching 
effect due to the flow field is neglected, which is more feasible for small molecules [21]. From the 
mathematical point of view, it brings great convenience since a weak maximum principle for |d| 
holds (cf. [221 Theorem 3.1], also [21]). When the stretching effect is taken into account, several 
results have been obtained. For instance, by taking 

Ai = -1, A 2 = 2a-1, 

Hi = 0, fi2 = -a, fj, 3 = l- a, m = v, 

fiFj = a(2a — 1), fj,Q = (a — l)(2a — 1), 

for some constants a G [0,1] and v > (cf. |4H Remark 3.2]), we can formulate a simplified 
version of the E-L system (|l.ip - (|1.3p such that the molecule has a general ellipsoid shape (cf. 
e -g-> [23U33] ). The parameter a is related to the shape of the liquid crystal molecules. For 
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instance, the spherical, rod-like and disc-like liquid crystal molecules correspond to the cases 
a = 2, 1 and 0, respectively (cf. [T4"l l33]). Under the above choice of coefficients, if a ^ \ then 
A2 7^ and the stretching effect is kept into account. We refer to [33] for the case with rod-like 
molecule (a = 1) , where existence of solutions in the regularity class 

uG L co (0,T;H 1 )r\L 2 {0,T;H 2 ), d G L°°(0, T; H 2 ) n L 2 (0, T; H 3 ) (1.10) 

was obtained in 2D, or in 3D under the assumption that the viscosity coefficient v is sufficiently 
large with respect to proper Sobolev norms of the initial data. See also [13] for certain blow-up 
criteria of the local classical solutions. For the general case a G [0, 1] , we refer to |HU28t l4"U] for 
detailed analysis on well-posedness as well as long-time behavior of the system. As was pointed 
out in [33] , due to the lack of maximum principle for the director equation when A2 7^ , it seems 
to be difficult to define (energy bounded) weak solutions in the usual sense as in [2TJ[22] and 
one has to deal with more regular solutions like (|1.1U|) . Recently, existence of certain suitably 
defined weak solutions with energy bound was obtain in [2], with the help of an appropriate 
choice of the test functions that leads to a rigorous weak formulation of the problem. This 
method has also been used in [9], where the authors prove the existence of weak solutions for 
a non-isothermal variant of the liquid crystal system in [2,23,28,33,40] with Neumann (for d) 
and complete slip (for u) boundary conditions. 

We note that the liquid crystal system studied in [21 EHl E31 SQ] is much simpler than the 
general E-L system (|l.ip - (|1.3p . because the six independent Leslie coefficients fj,±, [x% are 
only characterized by two parameters a, v . However, it has been shown that relations between 
those coefficients are crucial to retain certain basic properties of the general E-L system (|1.1|) - 
(|1.3p (with penalty), e.g., the dissipative basic energy law [22l|5l]. In particular, when the 
stretching effect is taken into account (i.e., A2 is not necessary to be zero), in [41] the authors 
provided sufficient conditions on the Leslie coefficients to ensure that the total energy of the E-L 
system (jl.ip - (jl.3p is dissipated (cf. Section 2.1) and they discussed the important role of the 
so-called Parodi's relation (cf. (|2.5p ) in the well-posedness and stability of the system. Besides, 
well-posedness and long-time behavior of E-L system (|l.ip - (|1.3p were obtained under suitable 
assumptions on the Leslie coefficients (for instance, the fluid viscosity /X4 is sufficiently large). 
They work within the same regularity class (cf. (jl.lOp ) for the solution (u, d) as in [li y 33 1 [40]. 
due to the same technical difficulty, i.e, the lack of control on the L 00 -norm of the director 
field d . It is worth mentioning that in the recent work [37] a sufficient and necessary condition 
on the Leslie coefficients was shown to ensure that the energy of the original E-L system with 
constraint |d| = 1 is dissipated, and local well-posedness as well as global well-posedness for 
small initial data were obtained. We refer to the recent work [3"lll2 t [T3" l [20, 36, 321S2] and the 
references cited therein for mathematical results on various (simplified) liquid crystal systems 
under the constraint |d| = 1. 

The first aim of this paper is to show that, when the stretching effect is taken into account 
(namely, A2 7= is allowed) and without any restriction on the size of the fluid viscosity as 
well as the initial data (as it is for the uncoupled 3D Navier-Stokes system), it is possible to 
obtain the existence of certain suitably defined weak solutions with finite energy for the E-L 
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system (jl.ip - (jl.3p (see Theorem 13. 3p . The key point relies on an appropriate choice of the 
test functions leading to a rigorous weak formulation of the system (cf . Definition 13. lj) . This 
technique is analogous to the one first used in [2] for the simplified system, but with more 
delicated estimates. In order to treat the nonlinear stretching terms as well as higher-order 
stress terms in the system, we make use of the a double-level approximate scheme that is a 
combination of a suitable Faedo-Galerkin approximation and a regularization procedure. In 
particular, a nonstandard but physically meaningful regularization of the momentum equation 
is introduced by adding to it a r-Laplacian operator acting on the velocity u, i.e., we add in 
the stress tensor a term of the type |Vu| r_2 Vu. Then we shall pass to the limit after obtaining 
suitable uniform estimates only from the basic energy law. 

In the second part of the paper, we will consider local classical solutions to the E-L system 
(|l.ip - (|1.3p . Under the assumption that the Parodi's relation (cf. (|2.5p ) is satisfied (namely, 
under the set of assumptions Case 1 below), we first prove the existence and uniqueness of 
local classical solutions to the E-L system (|l.ip - (|1.3p subject to periodic boundary conditions 
(cf. Theorem 14. 3p . The proof also indicates that the Parodi's relation plays an essential role in 
obtaining proper higher-order energy inequalities for the E-L system. Next, we address a Beale- 
Kato-Majda type criterion that characterizes the first finite singular time of the local classical 
solution in terms of ||V x u||£oo as well as ||Vd||^oo (cf. Theorem 14. 4p . For the uncoupled 
incompressible Navier-Stokes equations (or Euler equation), in their pioneering work pQ Beale, 
Kato and Majda showed the criterion in terms of L 00 -norm of the vorticity V x u. As far as 
liquid crystal systems are concerned, the situation is more involved due to interactions between 
the fluids and molecules. In the recent work [37] the authors obtain a BKM type criterion for 
the full E-L system with constraint |d| = 1 in terms of ||V x u||^oo and ||Vd||ioo (see also [12] 
for a simplified system). For our system (jl.ip - (jl.3p the main difficulty comes from the lack of 
control on the term ||d||£°° , which brings troubles to estimate the nonlinear kinematic transport 
terms as well as the higher-order Leslie stress tensors. 

Remark 1.1. With minor modifications of the proof of Section 5, we can obtain existence 
of global weak solutions also for the case of the homogeneous Neumann boundary condition 
for the director field, i.e., d n d = 0, which is suitable for the implementation of a numerical 
scheme (cf. [25]). The case of periodic boundary conditions for u and d can be treated with 
even simpler computations. However, when we consider classical solutions (and their blow- 
up criterion), we are only able to handle periodic boundary conditions. The main technical 
reason is that the periodic boundary conditions allow us to integrate by parts to obtain proper 
higher-order differential energy inequalities (cf. also [33][^0] .) 

Plan of the paper. The remaining part of this paper is organized as follows. Section 2 
contains the notations and the constraints on the Leslie coefficients. In Section 3, we first 
introduce a proper weak formulation of the E-L system and state the weak solutions existence 
theorem. After obtaining some useful (lower-order) a priori estimates from the basic energy 
law, we prove the existence of weak solutions with finite energy. Finally, Section 4 is devoted to 
the proof of local well-posedness of local classical solutions and a BKM type criterion. 
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2 Preliminaries and assumptions 



Let fi C M 3 be a bounded domain with smooth boundary T. We denote by L p (£l;R) and 
\Y m >P(p,-M.) the usual L p -spaces and Sobolev spaces of real measurable functions on fi and 
by L P (Q;R 3 ), L p (fi;M 3x3 ), W m ' p (Q; R 3 ) , W m ' p (n;R 3x3 ) the corresponding spaces of vector 
functions. Sometimes they will be simply denoted by L p , W m ' p . If / is an interval of R + and 
X a Banach space, we also use the function space L P (I;X), 1 < p < +oo, which consists of 
p-integrable functions with values in X . Next, we introduce some function spaces related to 
the Dirichlet boundary value problem (cf. |2iy34|): 

v = c^(n-,R 3 ) n {v : v • v = o}, 

H = the closure of V in L 2 (Q; R 3 ), 
V= the closure of V in Wq' 2 (Q; IR 3 ), 
V' = the dual of V. 



When we consider the periodic boundary value problem, we set fi = Q , where Q is a unit box 
in R 3 . Then we denote by H™(Q) , m G N, the space of functions which are in W^ 2 (R 3 ) and 
periodic in space with period Q. For the sake of simplicity, we denote the inner product on 
L 2 -spaces by (•, •) and the associated norm by || • || . Throughout this paper, the same letter C 
stands for a constant which may be different each time it appears. 

We consider suitable conditions on the physical coefficients Ai, A2, (jli, ■■■,(J>6 1 so that the E-L 
system (|l.ip - (|1.3p . together with proper boundary conditions, obeys certain dissipative energy 
law (cf. |22j). More precisely, we take the constraints 

Ai < 0, (2.1) 

/xi > 0, fi 4 > 0, (2.2) 

M5 + M6>0, (2.3) 

Ai = //2-A l 3, A 2 = /U5-/i 6 , (2.4) 

+ ^3 = A*6 - A*5- (2.5) 

Assumptions ()2.ip ~ (|2.2p are made to provide necessary conditions for the dissipation of the 
system (cf. [8l fl8t[T9] ) . Relations (|2,4p are necessary conditions to satisfy the equation of motion 
identically (cf. [181 Section 6]). 

Equation (|2.5p is called Parodi's relation, which is derived from Onsager reciprocal relations 
expressing the equality between flows and forces in thermodynamic systems out of equilibrium 
(cf. [29]). The Parodi's relation yields a constraint on the Leslie coefficients such that the dynam- 
ics of an incompressible nematic liquid crystal flow now involves five independent coefficients in 
(|1.5p instead of six. In [41] . the authors show that the Parodi's relation also serves as a stability 
condition for the E-L system (|l.ip - (|1.3p . 

In this paper, we will assume the following two different sets of hypotheses on the coefficients, 
as proposed in [H] . 
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Case 1 (with Parodi's relation). Suppose that (|2.ip - (|2.5p are satisfied. Moreover, we 
assume 

^-</i 5 + ^ 6 . (2.6) 
— Ai 

Case 2 (without Parodi's relation). Suppose that (j2. lj) — ()2.4[) are satisfied. Moreover, we 
assume 



|A 2 - M2 - M3| < 2^-X lx /fj, 5 + /i 6 . (2.7) 

As far as the potential is concerned, for the sake of simplicity, we will always assume that 
it satisfies the Ginzburg-Landau approximation, that is 

W(d) = ^(\d\ 2 -1) 2 . (2.8) 

Remark 2.1. In general, W may be written as a sum of a convex part and a smooth, but 
possibly non-convex one, for instance, 

W e C 2 {R 3 ), W>0, 

W = W 1 + W 2 s.t. Wi is convex and W 2 G C 1 ^ 3 ), VW 2 G C 0,1 (]R 3 ; R 3 ). 

With minor modifications in the proof as in [5] , it is easy to check that our result on existence of 
weak solutions to the E-L system (|l.ip - (|1.3p (cf. Theorem 13.31 below) also holds for the general 
potential W given above. However, in order to prove existence of (local) classical solution (cf. 
Theorem 14. 3p . additional assumptions on smoothness of the potential W should be added. 

3 Existence of weak solutions 

First we state the weak formulation of the problem under consideration, in which the momen- 
tum equation (jl.ip together with the incompressibility condition (|1 .2|) are replaced by integral 
identities. 

Definition 3.1. (Weak solution). A pair (u, d) with u(0, •) = Uo, d(0, ■) = do, a.e. in £1, is 
called a weak solution to the E-L system (jl.ip — (jl.3j) subject to the Dirichlet boundary conditions 
(|1.8p and (|1.9|) , if it belongs to the class 

u G L°°(0,r;L 2 (n;R 3 )) nL 2 (0,T; W 1,3 («;R 3 )), (3.1) 

dtu G L 2 (0, T; W' 1 '^ (0; M 3 )), (3.2) 
d G L oo (0,r ; ^ 1 ' 2 (O;lR 3 ))nL 2 (0,r;Ty 2 ' 2 (f];M 3 ))nif 1 (0,r;Lt(O;lR 3 )), (3.3) 
and, moreover, 

/ u(t, •) ■ V<p = 0, for a.e. t G (0, T), (3.4) 
Jn 

for any if G W l > 2 (n;R), 



(dtu, v) ,6 i 6 — / u Cg) u : Vvdx 
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/ (Vd Vd) : Vvdx - / a : Vvdx, for a.e. t G (0, T), 
Jn Jn 



(3.5) 



for any v G Wq' 6 (0; R 3 ) such that V • v = 0. Finally, the equation for the molecule director is 
satisfied in the strong sense, that is 



d t + (u • V)d -wd+ ^Ad 



1 

do|r 



(Ad - V d W(d)) , a.e. in (0, T) x 0, (3.6) 
a.e. on (0,T) x T. (3.7) 



Remark 3.2. A weak solution (u, d) corresponding to Neumann or periodic boundary condi- 
tions can be defined in a similar way. 

Introduce now the total energy of the system (|l.ip - (jl.3p , consisting of kinetic and potential 
energies, given by 

£ = l\M 2 + Jl|Vd|| 2 + [ W{d)dx. (3.8) 
2 2 Jn 

The main result of this section read as follows and its proof of will be postponed to Subsection 
3.2. 

Theorem 3.3. (Existence of weak solution with finite energy). Let C M 3 be a bounded 
domain of class C 2 . Assume (12.81) and the initial data such that 



u G L 2 (n;R- 
d G W 1,2 (£l; 



divuo = inL 2 (n), 
, do|r £ w'^' 2 (T; R' 



(3.9) 
(3.10) 



Then, for both Case 1 and Case 2, the Dirichlet problem i3.4\) - i3. 7j) possesses a global-in-time 
weak solution (u, d), in the sense of Definition 13.11 Moreover, for a.e. t G (0,T), the following 
energy inequality holds, namely, for Case 1, 



/'I"/" U|d T Ad| 2 + ^|Vu| 2 )(ix-^||Ad-V d W(d)|p 



r/r 



+ 



whiie for Case 2 



1^5 + ^6 



(A2) 2 



|^d|| 2 dr < 5(0), 



£(<) + f\l (/,i|d T ^d| 2 + ^|Vu| 2 ) dx + n(\\Ad\\ 2 + 

J _ J ^2 



dr < 5(0). 



3.1 A priori estimates 



We establish here a number of formal a priori estimates. These will assume a rigorous character 
in the framework of the approximation scheme presented in Subsection 3.2 below. 

Let (u, d) be a solution to the E-L system (|l.l|) - (jl.3|) subject to one of the following boundary 
conditions: (i) u|r = 0, d|r = do; (ii) u|r = 0, <9 n d|r = 0; or (iii) f2 = Q with u, d periodic 
in space. Following the detailed calculations as in |41| Lemma 2.1, Lemma 2.2] and |22l Theorem 
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2.1], for all the three types of boundary conditions for d (i.e., Dirichlet, Neumann and periodic), 
we are able to obtain some formal energy inequality for the E-L system (jl.ip - (jl.3p . More 
precisely, we have 

j£(t) = - [ (fn\d T Ad\ 2 + ^\Vu\ 2 + (^ + ^)\Ad\ 2 )dx 

+Ai \\Mf + (A 2 - /xa - /i 3 )(AA, Ad). (3.11) 

As a consequence, 

• Case 1. The total energy £(t) is decreasing in time and it holds 

p(t) = -^( m |d^d| 2 + ^|Vu| 2 )^ + i-||Ad-V d ^(d)|| 2 

-(^5 + Ai 6 + ^)pd|| 2 <0. (3.12) 

• Case 2. The total energy £ (i) is decreasing in time and there exists a small constant 
rj > such that 

j t £(t) = - J^l^Adf + ^lVu^ + i^ + ^Ad^dx 
+A 1 ||AA|| 2 + (A 2 -^ 2 -/x 3 )(AA,^d) 

< -J (MxId^dp + ^lVul^dx-^Pdf + IIATH 2 ) 

< 0. (3.13) 

Remark 3.4. We can easily see from the equation (|1.3p that the energy dissipations in Case 
1 and Case 2 are indeed equivalent. In particular, recalling also (|1.4p . we have 

||Ad|| 2 + ||Ad - V A W{d)f « ||Ad|| 2 + ||AA|| 2 . 

From now on we just present the case in which (u, d) satisfies the boundary conditions of 
Dirichlet type. Corresponding results for Neumann or periodic boundary conditions can be 
obtained in a similar way. 

From the previous energy laws we can obtain now (formal) a priori energy estimates for the 



solutions to the E-L system 

ll u llL°°(0,T;L 2 (n;R 3 )) < C, 1 1 " 1 1 ^2 ( ,X; Wq 1 ' 2 (fi;» 3 )) - C " ( 3 ' 14 ) 

1^11^(0^.^1,2(0^3)) < C, (3.15) 

|| - Ad + V d W(d)\\ L 2 {m 2 {n , R3)) < C, (3.16) 

ll^d|| L 2 (0iT;L 2(o ; R3)) < C, \\N\\l 2 (0,T;L 2 {Q;R 3 )) < C, (3.17) 

||d T Ad||x,2( 0) T;L2(n;R 3 )) < C i ( 3 - 18 ) 



where C is a positive constant depending on ||uo||, ||do||wi,2, 0, and the coefficients of the 
system, but it is independent of T. 
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Hence, estimates (|3.14p - (|3.18p imply that 

u G L co (0,T;L 2 (n;R 3 ))r\L 2 (0,T;W^ 2 (n;R 3 )), (3.19) 

d G ^(O.T;^ 1 ' 2 ^;! 3 )), (3.20) 

-Ad + V d VF(d) G L 2 ((0,T) xQ-R 3 ), (3.21) 

Ad, A/" G L 2 (0,T;L 2 (O;M 3 )), (3.22) 

d T Ad G L 2 (0,T;L 2 (S1;]R)). (3.23) 



It follows from (|3TT9l) that u <£> u G L 2 (0, T; L2 (fi ; R 3x3 )) . Then, using the incompressible 
condition (|1.2p . we get 

(u- V)u = V- (u®u) G L 2 (0, T; W -1 'i (O; R 3 )). (3.24) 
From the definition of W and (|3.2(jp . we have 

VdV^(d) G L 2 (0, T; L 2 (fi; M 3 )). (3.25) 
This together with (|3.2ip and the Dirichlet boundary condition for d yields that 

d G L 2 (0, T; W 2 ' 2 {£1; R 3 )), (3.26) 
where the bound also depends on Hdoll^ji 2 ^ r • As a consequence, we have 

(u-V)d, tod G L 2 (0,T;L§(O;lR 3 )). 

Then from (|3.22[) and the definition of M we deduce that 

d t d G L 2 (0, T; L§ (ft; M 3 )). (3.27) 

We recall the well-known anisotropic Sobolev embedding theorem for the Banach space 

V 2 ((0, r)xfl) = L°°(0, T; L 2 (ft)) n L 2 (0, T; W 1 ' 2 ^)) 

such that, when the spatial dimension is three, it holds 1^((0,T) x Q) <—$■ L~r((0,T) x Q) . 
Therefore, we infer from (|3.20p and ()3.26p that 

VdGL^(0,T;Lf(O;R 3x3 )). (3.28) 

In a similar manner, we get 

u G L% (0,T;Lir(O;E 3 )). (3.29) 

Besides, using the fact that for s G [0, 1] (see e.g., [24^ Definition 1.1, pp. 27] for the definition 
of interpolation space (•, -) s ), 

(L 2 (0, T; W 2 ' 2 {n; R 3 )), L°°(0, T; W 1,2 (fi; K 3 ))) s = L 1 ^ (0, T; VF 2 " s ' 2 (0; M 3 )), 
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and the Sobolev embedding in 3D such that W 2 S ' 2 (Q;R 3 ) ^ L^- 1 (O; R 3 ) , then, taking 
s = | , we have (cf. [30] ) 

(L 2 (0,T;W 2 ' 2 (fi;M 3 )),L oo (0,T;VF 1 ' 2 (O;]R 3 ))) s L 10 (0, T; L 10 (J7; M 3 )). 

As a result, we infer from (|3.20p and (|3.26p that 

d G L 10 (0, T;L 10 (S1;]R 3 )). (3.30) 

From the above estimates (|3.17p . (|3.18p . (|3.28p . (|3.30p and proper interpolation inequalities, we 
see that the nonlinear stress terms fulfil 

Vd© Vd, Ad® d, N®d G L 2 (0,T;Lt(ft;M 3x3 ))nLi((0,T) x Q;R 3x3 ), (3.31) 
(d T Ad)d®d G L 2 (0,T;Lt(O;M 3x3 ))nL^((0,T) x ^;M 3x3 ). (3.32) 

For instance, for the nonlinear stress tensor of the highest order (cf. (|3.32p ). we have used the 
following facts 

||(d r Ad)d ® d|| i2(oT;i 6 (n . R3x3)) < Hd^dH^^.^^.K))!^!! 2 ^^.^^.^)) 

< C\\d T Ad\\ L 2( 0tT . L 2(Q. R ))\\d\\ 2 LOO (q^ t . w i,2(q,- 

and (cf. (133U]) ) 

||(d T Ad)d®d|| L K, ((0jT)xC;R3><3) < ||d T >ld|| L 2 (0iT . L2(n . R)) ||d|| 2 1 o (0)T . il O (n: 

Finally, since 

L 2 (0,T;Lt(O;M 3x3 )) C L 2 (0, T; iS (fi; M 3x3 )), 
then the (distributional) divergence of the Leslie stress tensor <r satisfies 

V • cr G L 2 (0,T;lF- 1 't(O;lR 3 )). 

3.2 Proof of Theorem 

The proof of Theorem 13.31 consists of several steps. As in [2j[9], we shall construct a suitable 
family of approximate problems whose solutions weakly converge (always up to subsequences) 
to certain limit functions that solve the problem in the sense of Definition 13.11 

Step 1. Approximation. We introduce a double-approximation scheme that consists of 
a standard Faedo-Galerkin method with an approximation of the convective term as well as a 
regularization of the momentum equation (|3.4p - (|3.5p by adding an r-Laplacian operator acting 
on the velocity u . 

We take the orthonormal basis {v n }^ =1 of the Hilbert space V. Fixing M, N G N such 
that M < N, we consider the finite-dimensional space Xjy = span{v n }^ =1 . The approximate 
velocity field utv,j\/ G C 1 ([0, T]; Xn) solves the following equations 

-57 / ujv,M-vcfc- f [u/v,m]m ® ujv )M : Vvda; + \- f \Vu NtM \ r ~ 2 Vu NtM ■ Vvdx 
at Jn Jci M Jn 
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/ 

J CI 



Vd NtM Vd NM : Vvdx - / a N)M : Vvdx, for all t £ (0,T), (3.33) 

Jn 

/ uw,m(0, •) ■ vdx = / uo-vdx, (3.34) 
Jn Jn 

for any v £ Xjv and certain fixed r G (-jP, +oo) . Here the symbol [v] A / denotes the orthogonal 
projection onto the finite-dimensional space Xm = span{v n }*£ 1 . The r-Laplacian regular- 
ization term jj\VuN,M\ r ~ 2 ^^N,M in (|3.33|) is introduced to obtain enough regularity for the 
velocity field, which enables us to deduce certain energy inequalities for the limit functions 
(um^m) after we pass to the limit N — > +oo (see Step 2 below). The approximate function 
d-N,M for the molecule director is determined in terms of Uat,m as the unique solution of the 
parabolic system 

A 2 
M 

= Ad NtM - V d W(d N)M ), in (0, T)xQ, (3.35) 

d N ,M = d |r, on(0,T)xr, (3.36) 
djv.Af(0,-) = d o, infi, (3.37) 



where 



^Jv.Af = -(Vuat iA / + V t uat i m), un,m = -(Vuat iM — V t ujv,jw) 



For any fixed M and N, we can solve problem (j3.33[ >— (j3.37p by means of a fixed-point 
argument as in \10\ Chapter 3] (cf. also |21j for a simplified model of the E-L system). Indeed, 
we observe that all the a priori bounds derived formally in the previous section still hold for the 
approximate problem. Hence, if we fix u £ C([0,T];Xn) , then we can find d = d[u] solving 
(|3.35p - (|3.37|) . Inserting d[u] in system (|3.33p - (|3.34|) . we can define a mapping u h-» T[u] , where 
u = 7"[u] is the solution of the system. On account of the a priori bounds, we can easily show 
that T admits a fixed point by means of the classical Schauder's argument on (0, To), with 
< Tq < T . Finally, since the a priori estimates are independent of the final time To , we are 
allowed to conclude that the approximate solutions can be extended to the whole time interval 
[0, T] . The details are omitted here. 

Step 2: Passage to the limit as N — > +oo. Hereafter we just treat Case 1 since Case 
2 can be handled in a similar way (cf. Remark 13. 4p . From (|3.33|) - ()3.34p . we deduce that the 
regular approximate solutions (ujv,Af> djv,Af) satisfy 

h\u NM (t)\\ 2 + hvd NiM (t)\\ 2 + I w(d N>M (t))dx 

1 1 Jn 

+ J J (ml d T A N ^ M d N ,M | 2 + ylVu^A/ 1 2 + -^r|Vu A r jA /| r ^ dxdr 



[ \\Ad NjM -V d W(d NjM )\\ 2 dT 
M Jo 
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+ 



v 2 

+ + 



(V 2 



\AN.Md-N.M\\ 2 dr 



= ^l|uiv,A/(0)|| 2 + ^|Vd J v iA/ (0)|| 2 + f W(d NM (0))dx (3.38) 
1 1 Jn 

and from the formal estimates in the previous section that the following convergent results hold 
(up to a subsequence) 

u N:M ^u M weakly-(*) in L°°(0, T; L 2 (0; M 3 )) n L 2 (0, T; W 1,2 (fl; M 3 )), (3.39) 

d N ,M -> d M weakly-(*) in L°°(0, T; W 1 ' 2 ^; M 3 )) n L 2 (0, T; iy 2 ' 2 (ft; M 3 )), (3.40) 

d t d NiM -> <9 t d M weakly in L 2 (0,T; L§ (ft;M 3 )). (3.41) 

By virtue of (|3.40p . (|3.4ip and the Aubin-Lions lemma (cf. e.g., [32]), we have, for any small 

0<£« 1, 

d NM -> d A / strongly in C([0,T]; VF 1_ ^ 2 (0; M 3 )) n L 2 (0,T; VF 2 ~^' 2 (0; M 3 )). (3.42) 
Then a simple interpolation argument yields that 

d N)M -»• d M strongly in ^((0, T) x O; M 3 ), for 77 G (1, 10) , (3.43) 

Vdiv.ii/ -»• Vd A/ strongly in ^((O.T) x 0;]R 3x3 ), for 771 G (1, ^ J . (3.44) 



3 r 

On account of the regularizing r-Laplacian term introduced in (|3.33|) . we obtain the following 
additional regularity from the corresponding energy estimate 

^l|Vu J v,M||ir ((0j T)xQ;R3x3 ) < C, for r G (^y,+oo^ , (3.45) 

where the positive constant C does not depends on M and N . Hence, 

1 t r ( 10\ 

^IVu^mI^Vu^a/ G L—((0,T) x ft;M 3x3 ), with — - G ( . (3.46) 

Estimate (|3.45p implies that 

Vu7v,a/^Vu m weakly in L r (0, T; L r (ft; M 3x3 )). (3.47) 
Furthermore, we have (recall also (|3.31|) - (|3.32|) ) 

d t u NtM -> d t u M weakly in (0, T; VF" 1 '^! (0; M 3 )) + L 2 (0, T; W' 1 ' 2 ^; K 3 )). (3.48) 
Combining (|3.47|) with (|3.48p and using the Aubin-Lions lemma once more, we get 

ujv,m -> um strongly in L r (0, T; VF 1 ^' r (fi; M 3 )), < f « 1, (3.49) 
which together with the Sobolev embedding yields (recall that r > ^ > 3) 

U7v,Af -»• u M strongly in L r (0, T; L°°(0; M 3 )). (3.50) 
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Then interpolating (|3.39p with (|3.50p . we obtain 

u NM -> u M strongly in L r+2 ((0,T) x ft;M 3 ). (3.51) 
A combination of (|3.5ip with ()3.44p gives 

u N M ■ Vd N M -> u M • Vd M strongly in L Sl ((0,T) x 0;M 3 ), with — = — h — . (3.52) 

si r + 2 r?i 

Let r = 4^ + ? > and take r\\ = ^ — for certain small < < min{?, 1} . It is easy to 
see that 

^(1±^>^(A + A<),> 2 , 

provided that < ?' << 1 . Similarly, combining (|3.47p and ()3.43p . one obtains, for some s 2 > 2, 
A n ,m&n,M ->■ ^Afdjw, ojn,M^-N,M ->■ u M &M, weakly in L S2 ((0,T) x S7;M 3 ), (3.54) 
where 

The above relations (|3.52j) . ()3.54p and the L 2 -bound of — Adjv,Af + VdW^d^Af) (cf. estimate 
(133811 ) imply that 

d t d NtM -> d*d M weakly in L 2 (0, T; L 2 (ft; M 3 )). (3.55) 
Finally, we easily see from the strong convergence (|3.42p that (cf. assumptions on If) 

V d W(d NM ) -> V d W(d A/ ) strongly in L 2 ((0, T) x O; M 3 ). (3.56) 

As a consequence, we are able to pass to the limit as N — > oo in the approximate equation 
([335]) to get, a.e. in (0,T) x 0, 

A 2 

<9td M + UAf • Vd A/ - u M d M + —A M d M = Ad M - X7 d W(d M ), (3.57) 

M 

with 

d M = d , a.e. in (0,T) x T, (3.58) 
djvf(0, •) = do|r ; a.e. in Q. (3.59) 

In order to pass to the limit N —> +oo in the approximate momentum equation (|3.33p . we 
investigate the nonlinear terms. First, it follows from (|3.5ip that 

[utv,m]m Ujv,M -> [uju]Af u M strongly in L 1 ^ 2 ((0, T) x 0; ]R 3 ), (3.60) 

while (|3.44p yields, for 7? 2 £ (l, §) , 

Vd^A/ Vd NM -»• Vd M Vd A/ strongly in L^ 2 ((0,r) x 0;M 3x3 ). (3.61) 
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Next, from (fQ2|) and (pT51j) . we have 

^JV,Mdjv,M ® djv,M -> A M d M ® d M weakly in L 2 (0, T; L§ (0; M 3x3 )), 
wjv,Mdjv,M ® djv,M > -'A/d -i/ d \/ weakly in L 2 (0, T; L2 (0; R 3x3 )), 

while it follows from <^M), ([332]) and (|335|) that 

(ujvt,m • Vd^M) ® djv,M -> (um ■ Vd M ) ® d M strongly in L 2 (0, T; L§ (0; R 3x3 )), 
Stdjv.M ® djv,M <9 t d M ® d M weakly in L 2 (0, T; L^ 3 (Jl; R 3x3 )), m G M , | 

Concerning the stress term of the highest-order {6^ m An^m^-n,m)^n,m ® djv,M , we first infer 
from (jSMD and ilH that for some S3 G (1, 10/7) , 

d T NjM A N , M d N , M -> d^A M d M weakly in L S3 ((0,T) x n ; R). (3.62) 

This can be improved due to the L 2 -estimate of d^ m An,m&n,m (cf. (|3.38|) ) and uniqueness 
of the weak limit 

d^ )M ^JV,Mdjv,M -> d M A A /d M weakly in L 2 ((0, T) xfl;R), (3.63) 
Combining it with (|3.43j) again, we have, for some S4 G (1, 10/7) , 

{^n,m^n,m^n,m)^-n,m ® djv,M — >■ (d Af Ajv/d A /)dM ® dM, 
weakly in L S4 ((0,T) x ft;R 3x3 ). (3.64) 

Besides, estimate (|3.46[) implies 



|VuAr iM | r ~ 2 VuAr iM -> | Vuju | r_2 VuM, weakly in L~((0,T) x fi;R 3x3 ). 

Then the pair of limit functions {um,<^m) solves the problem (we note that in (|3.33|) the 
projection on Xm is kept in the convective term, see (|3.60p ) 

/ m m • Vpdx = 0, for a.e. t G (0,T), (3.65) 

for any <p G C°°(0) and 

/ (d t u M ,v)dT- [ [ [um] m ® um : Vvdxdr + — / / |Vu M | r " 2 Vu M : Vvdxdr 
Jo Jo Jn M Jo Jn 

= [ [ Vd M Vd M : V-vdxdr - [ [ <r M ■ Vvdxdr, for all t G (0, T), (3.66) 
Jo Jn Jo Jn 



for any v G C°°(0;]R d ) with V • v = 0, where 

&m = Mi (d M A M d M )d M ® d M + V-iNm ® dji# + // 3 d M ® Mm + ^A M 
+fi 5 (A M d M ) ® d M + M6dM ® (^Mdjif) 
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€ L^(0,T;L^(ft;M 3x3 )). 



(3.67) 



It remains to show that in (|3.66p 

\Vn M \ r ~ 2 Vu M = |Vu A /| r - 2 Vu A /. (3.68) 

On account of the U -regularity [r > ^) of Vu M (cf. (13171) ) and (pTBTD . we are allowed to take 
v = uj[j as a test function in (|3.66p . On the other hand, we multiply (|3.57p by Ad A / — VdW(djv/) 
and integrate on (0, t)xfl (this is possible since (|3.57p makes sense in L 2 (Qt) , cf. (|3.52|) - (|3.55p ). 
Adding the two resultants together, we get 

h\u M (t)\\ 2 + l\\Vd M (t)\\ 2 + [ W(d M (t))dx 
1 1 Jn 

+ J J (nM T M A M d M \ 2 + y |Vu A/ | 2 ) dxdr 

+ T7 I I |Vu A /| r " 2 Vu M : Vu M dxdr 
M Jo Jn 

+-!- / I ||Ad M -V A W{d M )\\ 2 dxdr 
A i Jo Jn 

i (\ ) 2 \ f t 

+ [H5+ + ^ \\A M d M \\ 2 dT 

= Jl|uAf(0)|| 2 + J||Vd M (0)|| 2 + I W(d M (0))dx. (3.69) 
z 1 Jn 

Passing to the limit as N — > +oo in (|3.38p and using the lower semi-continuity of norms, we 
infer from (|3.69|) that 

limsup / / \V\iN^M\ r dxdt < / |Vu A f | r ~ 2 Vu A / : VuMdxdt 
jv->+oo Jo Jn Jo Jn 

By means of the well-known Minty's trick for monotone operators (cf. [27]), we deduce that 
([3.68P holds. This concludes the passage to the limit as N — > +oo . 

Step 3: Passage to the limit as M — > +oo. The final step is to pass to the limit as 
M — > +oo in (|3.57p - (|3.59p and (j3.65p - (j3.67p . First, we observe that we are able to obtain the 
same convergence results as in (|3.39p ~ (|3.44p . (|3.56|) . (|3.61|) and (|3.64|) for u A / and d A / while 
letting M — > +oo, since the r-Laplacian regularization is not necessary to perform the limit in 
these terms. Moreover, we deduce from (|3.68p and (|3.69p that 

]gl|VUM||2r (( o,T)xn;M3x3 ) < C, r£ , (3.70) 

where C is independent of M . As a consequence, the following convergence result holds true 
(keeping in mind that -£j < ^ ) 

d t u M ^d t u weakly in L^t(0,T;VF _1 '^(O;IR 3 )). (3.71) 
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Using (|3.39p and the Aubin-Lions lemma again, we get 

u M -> u strongly in L 2 (0, T; PF 1_? ' 2 (ft; R 3 )). (3.72) 

Then we conclude that, for some S5 G (1,5/4) and S6 £ (1)5/3), 

u M • Vd M -> u • Vd strongly in L S5 ((0,T) x ft;M 3 ), (3.73) 

-4j\/dA/^^4d, w M djif ->• wd, weakly in L S6 ((0, T) x 0; E 3 ). (3.74) 
Next, from (|3.43p (for d.M instead of djv,jw) and (|3.74p . we deduce 

A M d M ®d M -> Ad®d weakly in L S7 ((0, T) xfi;M 3x3 ), 
w M d M ® d M ->■ ood (g) d weakly in L S7 ((0, T) x 0; M 3x3 ), 

with 57 G (1,10/7). Besides, for some ss G (1,10/9) and sg G (1,6/5), we have 

(uju • Vd M ) ® d M -> (u • Vd) d strongly in L S8 ((0, T) x 0; M 3x3 ), 
9td M «>dA/ -> d t d®d weakly in L 2 (0, T; L S9 (0; M 3x3 )). 

Finally, concerning the r-Laplacian regularization term, by the interpolation inequality and 
(j3.70p . we obtain 

M"^Vu M 

L r - 1 {(0,T)xn;R 3x3 ) 

< M"('- 1 )(''- 2 ) (M-r||Vu M ||L-((0,T)xn;K3x3)) l|Vu M ||iI ( O (0>r)><n . R3X 3 ) 

< CAfC'-iKr-a) with 0= , r ( r ~ 3 ) g (0,1), (3.75) 

(r — l)(r — 2) 

which yields 

M^Vu AJ ->• strongly in L r-1 ((0,r) x 0;M 3x3 ). (3.76) 

We are now in a position to pass to the limit as M — > +00 in (|3.57p ^ (j3.59p and (|3.65p - (|3.67|) . 
and finally recover the system (|3.4p - (|3.7p . The proof of Theorem 13.31 is complete. 

4 Local well-posedness and BKM type blow-up criterion 

In this section, we first prove the existence and uniqueness of local classical solutions to the E-L 
system (|1.1[) - (|1.3|) subject to periodic boundary conditions and then establish a BKM type blow- 
up criterion. For this purpose, the following results will be frequently used in the subsequent 
proofs 

Lemma 4.1. (cf. [IS])- F° r s > 0> there holds 

\\fg\\m < C(\\f\\ L ^\\ g \\ HS + \\f\\ H s\\ g \\ L ^), 

\\A s (fg) - fA s g\\ < C(\\f\\ H s\\g\\ L ~ + \\Vf\\ L ~ \\g\\ H .-i), 
where A is the Fourier multiplier such that Af(x) = X] fcgZ d(l + \k\ 2 )i e 2mk ' x f (k) . 
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Lemma 4.2. (cf. |35]). Let s > 0. Assume that T(-) is a smooth function on M. with T(0) = 
and f G H a (Q) n L°°(ft) . Tien we have 

\\F(f)\\ H s<c(l + \\f\\ L ^ +1 \\f\\ HS , 

where the constant c depends on sup ll-P 1 • 

0<fc<M+2, \y\<\\fh°° 

4.1 Local classical solutions 

Theorem 4.3. (Local classical solution). Let $7 := Q = (0, l) 3 be the unit cubic in M 3 . Suppose 
that the conditions in Case 1 are satisfied and s > 3 is an integer. For any uq G Hp(Q) with 
V • uo = and do G H p +1 (Q), there exists a Tq > depending on ||uo||#3, ||do||#4 and the 
coefficients of the system, but uniform in s, such that the E-L system (|l.ip - (|1.3p subject to 
periodic boundary conditions admits a unique local classical solution satisfying 

u G C(0,T ;H s p (Q)nL 2 (0,T ;H; +1 (Q)), (4.1) 
d G C(0, T ; H; +1 (Q) n T 2 (0, T ; H s p +2 {Q)). (4.2) 

Proof. In order to prove the existence and uniqueness of local classical solutions, we can first 
construct approximate solutions to the periodic value problem of the E-L system (jl.ip - (jl.3|) by 
the Galerkin method (e.g., [38] )• We note that here we no longer need the double-approximation 
scheme described in Section 3, since we can make use of higher-order estimates. After obtaining 
suitable uniform estimates on the approximate solutions, we are able to pass to the limit. Since 
the approximation procedure is standard, here below we just perform the necessary uniform 
higher-order energy estimates. 

Under the assumptions in Case 1, we can see that the (lower-order) uniform estimates 
(|3.14|) - (|3.18|) are still satisfied. Next, recalling |4H Lemma 5.1], under assumptions of Case 1 
the quantity A(t) = ||Vu|| 2 + ||Vd — VdW / (d)|| 2 satisfies the following differential inequality: 

j t A{t)<C{A{tf + A{t)), (4.3) 

where C is a constant that only depends on the coefficients fii, Aj, ||uo|| and ||do||#i . It follows 
from (|4.3p that there exists a certain T\ < +oo such that 

A(t)<C, tG[0,Tx], (4.4) 

for some constant C depending only on ||uo||#i, ||do||#2 and the coefficients of the system. As 
a result 

sup ||u(f)|| H i + ||d(t)|| H2 <C. (4.5) 

0<i<Ti 

In particular, by the Sobolev embedding theorem, we also have 

||d(t)|| L ~ <C, Vi€[0,Ti]. (4.6) 
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On account of estimates (|4.5[) and (|4.6p . we proceed to obtain higher-order estimates for u 
and d. For s > 3, applying A s to (jl.ip and testing the resultant by A s u, we get 

^l|A^|| 2 + f ||VA*u|| 2 
= -(A s (u • Vu),A s u) + (A s (Vd Vd), VA s u) 
-fi 1 (A s [(d T Ad)d d], VA s u) 
-fi 2 (A s (M d),VA s u) - /x 3 (A s (d M),VA s u) 
-fi 5 (A s [(Ad) d], VA s u) - fM 6 (A s [d (Ad)}, VA s u) 

i=l 

Moreover, applying VA S to ()1.3p and testing the resultant by VAM, we obtain 

^IIVA-df-i-HAAMf 

= -(VA s (u • Vd), VAM) - (A s (wd), AA s d) 

+ ^ (A s (Ad) , AA s d) + J-(A s [V d iy(d)], AA s d) 
Ai Ai 

11 

i=8 

For , using the incompressibility of u , we have 

K x = -(A s (u-Vu)-u-VA s u,A s u) 

< C||A s (u- Vu) -u- A s Vu||||A s u|| 

< C||Vu|| L oc||u||^ s . 

The term K<i can be simply estimated as 

K 2 < ||A s (Vd0Vd)||||VA s u|| 

< KllVA^f + CllVdlll^HVAMII 2 . 

Similarly to K\ , for K$ , we get 

Kg = -(VA s (u- Vd) -u- V(VA s d),VA s d) 

< C(||u|| H .+i HVdlli- + || Vu|| L oc ||Vd||^)l|VA s d|| 

< k||VA s u|| 2 + C(l + ||Vu|| L oc + ||Vd||| 00 )||VA s d|| 2 . 

By the interpolation inequality 

||Vu|| L «» <C||u|||r^||A s u||^, s>3, 
and since 2 3 _ 2 < | when s > 4, we have 

lldl&.IIAHic <q|d|| fl3 ||d|| H 4||Vu||l <JO ||u||| 1 ||A 3 u|| J s = 3, 
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l|d||^||A|||oo < C||d|| HS -i||d|| HS+ i||u|| Hl ||A s u||, a > 4. 
From (|4.5p we infer, for s > 3, 

lldll^PHloo < C(l + ||Vu||f 00 )(||A'u|| 2 + Hdll^lldll^O, t G [0,^]. (4.7) 

Then for the stress term A3 , noticing that is antisymmetric and using (|4.5p , (j4.7[) , we obtain 

K 3 = -/ii(A s [(d T Ad)d«)d],A s (A + w)) 
= -^i(A s [(d T Ad)d®d],A 3 ^) 

= -/ii||d T AMd|| 2 + /ii(A s [(d T yld)d® d] - [(d T A s id)d®d],A s i) 

< -^i||d T A s Ad|| 2 + K||VA s u|| 2 + C||A s [(d T Ad)d®d] - [(d T AMd)d® d]|| 2 

< -/ii||d T A s Ad|| 2 + K||VA s u|| 2 
+C(||d||£„||d||iH|A|| i oo + ||d|||o ||Vd|| I; o ||A|| ifs - 1 ) 2 

< -/xi||d T AMd|| 2 + k||VA s u|| 2 + C(l + ||Vu|| L oc)(||A s u|| 2 + ||d[|^ s _ a ||d||^+i) 
+C||Vd|| 2 oc ||A*u|| 2 

< -^i||d T A s Ad|| 2 + K||VA s u|| 2 

+C(1 + ||Vu|| L oc + ||Vd||| 00 )(||A s u|| 2 + HA'-MfllA^dll 2 ). 

Using the definition of A, u , and the fact that u) is antisymmetric, then it holds 

K± + K 5 + K 6 + K 7 
= -0*2 + /i 3 )(A s (A/ d), AM) - (fji 2 - fi 3 )(A s (M d), A s uj) 

-0*5 +/i 6 )(A s ((Ad) ®d),A s A) - 0t 5 -/i 6 )(A s ((Ad)»d),A^) 
:= P 1+ P 2+ P 3 + P 4 . 

On the other hand, using the equation (jl.3p such that \\Af + X 2 Ad = — Ad + VdW^(d) , we have 

Kg + K W 

= X 1 (A s M,A s {iod)) + X 2 {A s (Ad),A s {ud)) 

-A 2 (AW,A s (Ad)) - ^-(A s (Ad),A s (Ad)) 

-(A s (V d IU(d)), A>d)) + ^(A s (V d ^(d)), A s (Ad)) 
10 

:= 

i=5 

Recalling the condition (|2,4p and the Parodi's relation (|2.5p . we can recombine the above terms 
Pi,-.., Pio and eliminate the term with highest-order derivative as follows 



P1 + P7 + P 



10 



20 



= A 2 ( A s (TV ® d), A S A) - A 2 (AW, A s {Ad) ) + ^(A s (V d W(d)), A s (Ad) ) 



= A 2 (A s W ® d) - AW ® d, AM) - A 2 (AW, A s (Ad) - A s Ad) 



= A 2 (A S (A/" ® d) - AW ® d, A S A) + ^-(A s {Ad), A s (Ad) - A s Ad) 



+-^(A s Ad, A s (Ad) - A s Ad) + ^(A s V d ^(d),A s id) 
Ai Ai 

:= R\ + i? 2 + i?3 + i?4. 

p 2 + p 5 + p 9 

= Ai(AW,A s (wd)) - Ai(A s (A/"®d),A s w) 

-(A s (V d W(d)),A s M)) 
= Ai(AW,A s (wd) - A s wd) - Ai(A s (AA^d) - AW®d,A s w) 

-(A s (V d W(d)),A>d)) 
= -Ai(A s (AA ® d) - AW ® d, A s uj) - A 2 (A s (Ad), A s (wd) - A s wd) 

-(A s Ad, A s (ud) - A s ud) + (A s (V d ^(d)), A s wd) 

: = i?5 + i?6 + R7 + Rs • 

P 4 + P 6 = -A 2 (A s ((Ad)®d),A s u;) + A 2 (A s (Ad),A>d)) 



^3 + ^8 = -(fi 5 + n 6 )(A s ((Ad)®d),A s A)-^(A s (Ad),A s (Ad)) 



As a consequence, in order to estimate the terms K^, K§, Kq, K7, Kg, Kiq , we can turn to esti- 
mate the new terms R\, R13 . 
By the following bounds 



+^(A s (V d W(d)),A s (Ad)) 



-A 2 (A s ((Ad) ® d) - A s (Ad) ® d, A s w) 
+A 2 (A s (Ad), A s (ujd) - A s ud) 

Rg + Rio- 




-(Me + W;)(A s ((Ad) ® d) - A s (Ad) ® d, AM) 
+(Mfi + ^)(A s (Ad), A s (Ad) - AMd) 

R\l + i?i 2 + R13. 



d\\ H3 <c||d||| 4 ||vd||2 



1 1 



Ad|| L oo <C||d||| 3 ||d 
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we deduce, for t £ [0, Ti] , 

lldl&.HAdlliao < C||d||^ 3 ||d|| H 4 < C||Vd|||oo||d||^ 4 , when s = 3, (4.8) 
lldll^HAdllioo < C||d|| H .-i||d|| ff ^i||d||H3||d|| H 4 

< C||d||^_ 1 ||d||^ +1 , whens>4. (4.9) 

Next, using equation (|1.3[) . Lemmas 14.11 14.21 and estimates (|4.5p - (j4,9p . we have 

R x < K||AM|| 2 + C7(||d||^||AA||L- + ||Vd|| LO o||AA||^i) 2 

< K \\A 3 A\\ 2 + CUdlllr.CPIIioclldllioc + ||Ad||£„ + ||V d ^(d)|||cc) 
+C||Vd||| 00 (P||^_ 1 ||d||| 00 + IIAHioolldll^ + IIAdHlr..! + ||V d W(d)||| s _x) 

< k||A 3 A|| 2 + C(l + HVuHioo + ||Vd||| 00 )(||A s u|| 2 + ||d||| r .-i||d||^x) 
+C||d||^- 1 ||d||^ +1 (||Ad|| ifl ||Ad|| /f2 + l) 
+C||Vd||| 00 (||A s u|| 2 + HA^df + IIA^MH 4 + ||A s+1 d|| 2 ) 

< k\\A 3 A\\ 2 + C(l + ||Vu|| L oo + ||Vd|||o )[||A s u|| 2 + (1 + Hdll^-OUdll^+i]. 

R2 < C(P||^||d|| L oo + [|A|U«[|d||H.)(l|d||H'l|^lk- + l|Vd|| L ooPII^-0 

< k\\A s A\\ 2 + C||d[|^||Vu||ioc + C||Vd||£ c||A s u|| 2 

< k\\A s A\\ 2 + C(l + ||Vu|U« + ||Vd||| 00 )(||VA s d|| 2 + ||A s u|| 2 ). 

i? 3 < K||AA s d|| 2 + C , (||d||^P||L- + ||Vd|| Loo p|| Hs _ 1 ) 2 

< k\\A s A\\ 2 + C(l + ||Vu|| L « + ||Vd||| o)(||VA s d|| 2 + ||A s u|| 2 ). 

Ra < «||A^4|| 2 + q|A s V d W(d)|| 2 ||d|l!~ 
< k||A 3 A|| 2 + C||AM|| 2 , 

where for R^ we have used Lemma H~2l and (|4.6p . The term R§, Rq, R7, Rg can be estimated 
in the same way as R± , R2 , R3 , Ra , respectively. Besides, the term Rq,R\2 can be estimated as 
Ri while i?io,-Ri3 can be estimated as R2 ■ Moreover, from (I2.6P it follows 

Rn < 0. 

Finally, K\\ can be estimated similar to -R4 so that 

K n < K||AA s d|| 2 + C||A s d|| 2 . 
Collecting the above estimates together, taking k > small enough, we conclude 

^-(||A s u|| 2 + ||VA s d|| 2 ) + ^||VA s u|| 2 - -?-||AAM|| 2 
at 2 Ai 

< C(l + ||Vu|| L oo + HVdlllcc) (1 + ||A s u|| 2 + (1 + lldll^-OllVAMf) , (4.10) 
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for s > 3 and t G [0, Ti] . Here C is a constant depending on ||uo||#i, ||do||#2, T\ and the 
coefficients of the system. 
Denote 

y s (t) = ||A s u(i)H 2 + l|VA s d(t)|| 2 . 
When s = 3, using the Sobolev embedding theorem, from (|4.10p and (|4.5p we infer 

jY 3 (t) < C(l + ||A 3 u|| 2 + ||VA 3 d|| 2 )y 3 (t) < C(l + Y 3 (t) 2 ). (4.11) 

As a result, there exists a time To G (0,Ti] depending on y(0) (i.e., ||uo||#3 and ||do||#4) such 
that ^(i) < +oo. Moreover, we get 

sup \\u(t)\\ H3 + \\d(t)\\ H 4<C, 

0<t<T 

which, together with an application of the Sobolev embedding theorem, implies 

sup ||Vu|| LO o + ||Vd(t)|| L oo. (4.12) 

0<t<T 

Then we can use a simple induction argument. We have shown that Y^t) < +oo on [0, To]. 
Suppose, for k > 3, we have Yfc(i) < +oo on [0, To]. Then from (|4.10p and (|4.12|) we get 



J t Y k+i(t) < C(l + ||Vu|| L ~ + HVdllico) (l + ||A fc+1 u|| 2 + (1 + ||d||^)||VA A 



^ +1 u|| 2 + (l + ||d||^)||VA fc+1 d| 
< C(l + Y k+l (t)), (4.13) 

which yields 

Y k+1 (t) < (Y k+1 (t) + l)e c \ t G [0,T ]. 
Thus, for any integer s > 3, uo G H^{Q) with V • uo = and do G H^ +1 {Q) , we have 

sup ||u(t)|| H « + \\d{t)\\ Hs +i < C, 

0<t<T 

where C depends on ||uo||h s > 11*^0 II J e _s + 1 and the coefficients of the system, while the existence 
time To > only depends on ||uo||#3, ||do||#4. As a consequence, we obtain a priori esti- 
mates which are uniform for the approximate solutions on [0, To] . By a standard compactness 
argument, we obtain the existence of the local classical solution on [0,To] : 

uGT°°(0,To;iI p s (Q)nL 2 (0,To;F p s+1 (Q)), d G L°°(0, T ; ^ +1 (Q)nL 2 (0, T ; H s p +2 (Q)). (4.14) 

The uniqueness of the solution can be easily derived by using the Gronwall inequality in 
the norm in L 2 x H 1 , since we have gained enough higher-order estimates of the solution. 
Concerning the time continuity property of u, d, we can argue as in |38| . For any 5 > 0, we 
can take N = N(5) such that 

sup ^2 2 ^||A,u(t)|| 2 <^, 
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where Aj is the Littlewood-Paley decomposition operator (see, e.g., [38J Appendix]). Then for 
any t G (0, To) and <r such that t + <r G [0, To] , we have 

\\u(t + a)-u(t)f H . 

N 

< V 2 2 ^||A J u(t + a)-A j u(t)|| 2 + 2 sup V2 2 ^||A i u(t)|| 2 
; i ' ' * '-, v 

< |a| £ 2 s * / \\ dt u(t)\\ 2 dt + - 

j=-l J o 

< \a\(N + l)2 2NsN [ T °\\d t u(t)\\ 2 dt+ 5 -. 

Jo 1 

We easily infer from (|4.14p and the equation (jl.ip that u G T 2 (0, To; L 2 (Q)) . As a result, for a 
small enough, we have ||u(t + a) — u(t)|| 2 ^ 3 < 5 . The time continuity of d can be shown in a 
similar way. The proof is complete. I 

4.2 Blow-up criterion 

Theorem 4.4. (BKM type blow-up criterion). Suppose that the assumptions of Theorem 14.31 
are satisfied. Let (u, d) be the local classical solution to the E-L system ([l.ip — (|1.3p subject 
to periodic boundary conditions. Let T* be the maximal existence time of the solution. If 
T* < +oo , then 

[ (||Vxu(t)|| L » + ||Vd(t)||i 00 )<tt = +oo. (4.15) 
Jo 

In order to prove Theorem 14.41 we need the following lemma. 

Lemma 4.5. Suppose that the conditions in Case 1 are satisfied. For any uq G H^{Q) with 
V • uo = and do G Hp(Q) , M > and Tq > 0, let (u, d) be a local classical solution to the 
E-L system (|l.ip - (|1.3|) subject to periodic boundary conditions. If the following condition is 
satisfied: ^ 

I ° (||V x u(t)|| L .» + ||Vd(t)||ioo) dt < M, (4.16) 
Jo 

then 

sup (\\u(t)\\ Hl + \\d(t)\\ H2 )<C, (4.17) 

0<t<T 

where C is a constant depending on ||uo||#i , ||do||#2 , M and the coefficients of the system. 

Proof. In order to obtain the estimate (|4.17p . we shall derive some energy differential inequal- 
ity like (|4.10p . Observe that under assumptions in Case 1, uniform estimates (|3.14p — (|3.18p still 
hold. However, here the situation is more delicate because these estimates are rather weak and, 
in particular, we lose the control of ||d||^oo(o 5 T;L°°) • 

Then, applying the curl operator Vx to (|l.ip . we obtain 

(V x u) t + (u • V)V x u = (V x u) • Vu - V x (VdAd) + V x (V • cr), (4.18) 
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where we have used the well-known identity V • (Vd Vd) = V(|| Vd| 2 ) + VdAd and the fact 
that V x (V-) is the null operator. Besides, we recall the identity 

V x V x u = V(V • u) - Au, 

which together with the incompressibility condition (j 1 . 2 [) implies that 

V x V x u = -Au. 

Testing (|4.18p by V x u , we obtain that 

Il||Vxu|| 2 + ^||Au|| 2 

= ((V x u) • Vu,V x u) + (VdAd, Au) - ^(V • [(d T Ad)d ® d], Au) 
-/i 2 (V- (Af <8>d),Au) -M3(V • (d (g> AT), Au) 
-/i 5 (V • [(Ad) CS> d], Au) - /i 6 (V • [d ® (Ad)], Au) 

7 



m=l 

where we have used the fact 



((u • V)V x u, V x u) = - [ (u • V)| V x u\ 2 dx = 0. 

2 JQ 



Next, applying A to (jl.3p and testing the resultant by Ad, we get 

Ii|Ad|| 2 -^-||VAd|| 2 
2dt u 11 Ai" 11 

= -((Au • V)d, Ad) - 2 / ViUjViVjdkAdkdx 

JQ 

+(A(wd), Ad) - ^ (A (Ad), Ad) + ^-(AV d W(d), Ad) 
Ai Ai 

12 
m=8 

Obviously, we have the cancellation /2 + ^8 = . Since the Riesz operators are bounded in I? 
and Vu = (- A) _1 V(V x V x u) , we deduce || Vu|| < C\\V x u|| . Then for Ii , we see that 

h < l|V 

X U roc ||Vu||||V x u|| < C||V ||V x ull . 

Moreover, using the Gagliardo-Nirenberg inequality and the Holder inequality, for Ig and I\2 
we have 

h < C||Vu||||V 2 d|| 2 4 

< C||V x u||(||Vd|| L oc||VAd|| + ||Vd|||oo) 
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< K||VAd|| 2 + C||Vd|||oo||V xu|| 2 , 



Jia = - / V(|d| 2 d): VAddx + -||Ad|| 2 

£ JQ £ 

< C||d||i 6 ||Vd|| L 6||VAd|| + i||Ad|| 2 

< K||VAd|| 2 + C||d||J r i||V 2 d|| 2 + ^||Ad|| 2 

< K||VAd|| 2 + C||Ad|| 2 . 

Using integration by parts and the fact that u is antisymmetric, we deduce from the Holder 
inequality and Young inequality that 

I 3 = m ([(d T Ad)d®d],AVu) 
= /x x ([(d T Ad)d®d],A(A + a;)) 
= /ii([(d T ^d)d®d],A^l) 

= -Hi I V i(d k A kp d p didj)V lAijdx 



-Hi / (d k d p ViA kp ) 2 dx 
Jo 



~Hi / A kp (Vid k d p + d k V d p )didjV lAijdx 
JQ 

-Hi \ A kp d k d p {diV 'idj + V 'ididj)V ' lA^dx 
JQ 

< "TT / {dkdpViA^fdx + k||Au|| 2 

+C m (||d|| 2 00 + ||d||l 00 )||Vu|| 2 ||Vd|| 2 00 . 
Next, we infer from the following inequality 

||d|Uoo <C7||Vd||| 00 ||d||| 6 + C7||d|U6 



that 



|d|| 2 ^ + ||d||| 0O )||Vu|| 2 ||Vd|| 2 o 



c 



||Vxu|| 2 ||Vd|| 2 



l+givdii^Hdii^ + ciidib 

< C(l + || Vd||i<x,) II V x u|| 2 . 
As a consequence, we have 

h < ~J (4d P V^ fep ) 2 dx + K ||Au|| 2 + C/i 1 (l + ||Vd||l 0O )||Vxu|| 2 . 
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Concerning I4, , using integration by parts, we obtain (cf. e.g., [H] Appendix]) 

h + h = O2 + Ms) / djNiAAijdx + fa - Ms) (A", Aw d) (4.19) 
JQ 

and 

I 6 + / 7 = -(/U5 + /x 6 ) / IdjV/Ajil 2 ^ - (^ 5 + /i 6 ) / V idjd k A ki V lAijdx 

JQ JQ 

-(M5 + M6) / djV id k A ki V lAijdx + (// 5 - // 6 ) (^d, Awd) 
:= J1 + J2 + J3 + (M5 - Me) (-4d, Awd) . (4.20) 
On account of equation (|1.3p and integrating by parts, we get 

/10 = (Ad,Awd) + 2/ AdiViLUijV idjdx + (Ad,wAd) 

JQ 

= -Ai / djMiAujij dx - \ 2 {Ad, Awd) + (X7 d W(d), Awd) 
JQ 

+2 / Adi Viojij Vidjdx + (Ad, u Ad) 
J Q 

= - Ai (A/\ Aw d) - A 2 {Ad, Awd) - [ V / [(V d VF(d)) i d i ]V i ^ j dx 

JQ 

V [AdiUJijV 'idjdx + / AdjV^WjjV/djda; 
Jo 



and 



- Ai (A/", Aw d) - A 2 (Ad, Aivd) + J 4 + J 5 + J 6 , (4.21) 



r n = \ 2 (jV,A(Ad))+^-(Ad,A(Ad))-^(V d W(d),A(Ad))) 

= A 2 / MiAAijdjdx + 2A 2 / A/iV/AyV^cte + A 2 (A/", AAd) 
JQ ' JQ 

( A 2) 2 /" 1T7 / ^ J M2j„ , A 2 



|V i (^ i d i )| 2 ^ + T i (W d W(d),V(yld)) 

Al JQ A l 

A 2 / MiAA^djdx + A 2 / MiV lAijV idjdx - A 2 / V 'jAZ^A^-V 'idjdx 

JQ 

(A 2 j 



\2 /■ /\ ^2 



A l JQ A l JQ 



/ V ^ JJ ^A ife V z 4^ + ^ 2 - / Vi(V d W(d)) i V^yd i dx 

Al Jq Ai jq 

+^ / vKv d ^(d))^v^x 

Ai Jq 

A 2 / MiAA^dx - [ \djViAji\ 2 dx 

JQ Ai 
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A 2 f A 2 f 

+— / V 'lAdiAijV idjdx — — / AdjV lAijV idjdx 

x iJq x iJq 

f VlAijdjAikVldkdx + h f ^i[(y d w(d))idj]ViAijdx 

Ai Jq Ai 

:= A 2 / MiAAjjdjdx + J 7 + ... + Jn. (4.22) 
Jq 

By the Parodi's relation (|2.5p and the condition ()2.4p . the first term of the right-hand side of 
(|4.22p cancels with the first term of the right-hand side of (|4.19p . Besides, due to the relation 
()2.4p . the first term on the right-hand side of ()4.2ip cancels with the second term of the right- 
hand side of (|4.19p and the second term on the right-hand side of (|4.2ip cancels with the fourth 
term of the right-hand side of (|4.20|) . 

Keeping in mind these special cancelations between the nonlinear terms of the highest-order 
derivatives, we are able to estimate the remaining terms in ho, In , namely, Ji, J12 • It 

follows from (j2.3j) and (j2.6j) that J% + J7 < 0. Next, we have 

J2 + J3 < C(/U 5 + /i 6 )||Au||||Vu||||Vd||Loo ||d|| L oo 

|2 1 /~tf.._ 1 ,. . \2 / ||V7J II 3 IIJII3 1 IIJII \ IIV7JII2 IIV7^,,l|2 



< K||Au|| z + C(/i5 + / u 6 )^||Vd||2 00 ||d||2 6 + ||d|| L 6j ||Vd||£oc||V x u| 

< k||Au|| 2 + C(/i 5 + /ue) 2 (l + HVdllfoo) ||V x u|| 2 , 



J 4 < C||Au||||Vd|| L 6(||d|| L 3 + ||d||3 9 ) 
£ h\ Au|| 2 +c( 1 + HVdllloo ] ||Ad|| 2 . 



For J5, Jq, by an integration by parts and the fact that uj is anti-symmetric, we get 

J 5 + J 6 = 2/ AdiV lUJijV idjdx < k||Au|| 2 + C|| Vd^ || Ad|| 2 . 
Jq 

Integrating by parts and using the Holder inequality and the interpolation inequalities 

||Ad|| L 4 < CfllVAdll^lVdllloc + ||Vd|| L oo), 
||Vd||ioo < CUVdll^HVdllla, 

then we get 

A2 f 2A 2 f 

Js + J9 = ~T~ AdiAijAdjdx — / AdiV iA^V idjdx 

Ai Jq Ai J q 

< |A 2 |||Vu||||Ad|| 2 : 4 + |A 2 |||Au||||Vd|| L oo||Ad|| 

< k||Au|| 2 + K||VAd|| 2 + C(A 2 ) 2 ||Vd||| o(||V x u|| 2 + ||Ad|| 2 ) 
+C||Vd||| oc . 
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For Jio , we have 



(A 



2 



\2 



Jio < — r- — 1| Au|| ||d||x,oo || Vu|| || Vd||^oo 
Ai 

< k||Au|| 2 + C(A 2 ) 4 ||d||| 00 ||Vd|||oo||V x u|| 2 

< k||Au|| 2 + C(A 2 ) 4 (1 + ||Vd||| co )||Vd|||o ||V x u|| 2 . 

Finally, J\\ can be estimated just as J4 . 

Collecting the above estimates and taking k > small enough, from the the Young inequality 
and (|2.6p (i.e., | A2 1 can be bounded by ^5 + hq > 0), we infer 

j t (||V x u|| 2 + ||Ad|| 2 ) < G(t)(l + ||V x u|| 2 + ||Ad|| 2 ), (4.23) 



where 



G(t) = C 



1 + ||V x u|| Loo + ||Vd|||oo + (/i 5 + /i 6 ) 2 ||Vd||2 00 + )Ui||Vd||to 



(4.24) 



Here C is a constant depending on ||uo||, ||do||//i and the coefficients of the system. Finally, 
by the Gronwall inequality, we deduce 

|| V x u(t)|| 2 + ||Ad(i)|| 2 < (1 + ||V x u || 2 + ||Ad || 2 ) exp (J G(r)dr^) , t G [0,T ]. 

This and the lower order estimates (|3.14p - (|3.15p yield the conclusion (|4.17|) . The proof is 
complete. I 

Proof of Theorem 14.41 We prove the theorem by a contradiction argument. Suppose that 
(|4.15p is not true. Then there is a positive constant M such that 

/ (||Vxu(t)|U« + ||Vd(t)||ioo)£it<M, (4.25) 
Jo 



sup (||u(t)||^ + ||d(t)|| H2 )<C (4-26) 

0<t<T* 



which yields (cf. Lemma |4. 5 



where C is a constant depending on ||uo||#i , ||do||ij-2 , M and the coefficients of the system. 

Recalling the proof of Theorem 14.31 and using the bound (|4.26p . we are able to derive the 
following differential inequality (|4.10p for the quantity Y%(t) 

jY 3 (t) < C (1 + ||Vu|U=o + ||Vd|||oc) (Y 3 (t) + 1). (4.27) 
Combining the critical logarithmic Sobolev inequality (cf. pQ) 

||Vu||loo < C (1 + ||V x u|| + ||V x u\\l°° hi(e + ||u|| H 3)) , 
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with estimate (|4.26p . we get 



||Vu||l°° < C (1 + ||V x u||loo ln(e + ||u||#3)) . 
Then, from (|4.27p we infer 

j t {Y*{t) + e + 1) < C(l + ||V x u|| L oo + ||Vd|||ao)(y 3 (t) + e + 1) ln(F 3 (i) + e + 1). 
By using the Gronwall inequality, for t £ [0, T*), we conclude 

Y 3 (i) + e + 1 < exp fln[F 3 (0) + e + l]exp (c J (1 + ||V x u|| L oo + [|Vd||£,»)dT^ ^ . 

Therefore, if T* < +oo and (|4.25p holds, an application of the Young inequality gives 

I (1 + ||V x u\\ L oo + ||Vd|||oo)di < / (1 + ||V x u|| L °o + ||Vd||£cx>)cft < +oo. 
Jo Jo 

Then we deduce that Ys(t) is bounded for t G [0, T*) and the local classical solution (u, d) can 
be extended beyond t = T* . This leads to a contradiction of the definition of maximal existence 
time T* . The proof is complete. I 

Remark 4.6. Comparing our blow-up criterion (|4.15p to the one obtained in |12j for the sim- 
plified liquid crystal system with constraint |d| = 1 such that 

f (||Vxu(f)|| L o° + ||Vd(*)||£«»)dt = +oo, (4.28) 
J o 

they differ only on the power of || Vd||£°o . Observe that in [12] the term ||Vd||£°° as to be 
included in the criterion mainly because of the difficulty due to heat flow of harmonic maps. 
For the full E-L system with constraint |d| = 1 (cf. |37j). higher-order stress tensors will not 
introduce further troubles and the same criterion (|4.28p still works. However, in our case, the 
power of 1 1 Vd 1 1 ^oo has to be larger than those in [12^ 37]. due to the lack of control on [|d||nx> . 
Furthermore, the expression (|4.24p indicates that for different choices of the Leslie coefficients 
(always under assumptions in Case 1), the blow-up criterion (|4.15p can be improved. For 
instance, under the special case [i\ = /X5 + n§ = (e.g., taking a = ^ in the simplified liquid 
crystal system studied in ) , (|4.15|) can be replaced by (|4.28|) . In this case, we have also 

A2 = due to (|2.6p . namely, the stretching effect is neglected in the E-L system while the 
rotation effect is kept. When a = 1 , we refer to [13] for another type of criterion only in terms 
of the velocity 

T * Wu(t)\\ r Tp 2 3 

1 + ln(e + ||Vu(i)|| LJ >) r p 

but with a restricted range of p instead of the natural one | < p < 00 . 

If we neglect the parallel transport of the director d (i.e., the term — o;d + -j^-^d in (|1.3p ) 
and drop all the other Leslie stress terms in <r except the fluid viscosity term 114A, then the 
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E-L system f|l.ip f|1.3[) will be reduced to the simplest model studied in [21]. For this case, a 
logarithmical improved BKM's criterion in terms of the velocity field u can be obtained (cf. |26j), 
meaning 

IV x u(t)\\ BM0 =M < + ^ 



t - 

Jo i/M e + || V X u(l)\\BMO) 
We recall here the definition of BMO (Bounded Mean Oscillation) space 

BMO = J / € L /oc (M 3 ) : \\f\\ B MO = sup — / \f(y) - f BR{x) \dy < oo 
[ R>o,xm 3 \-t>R{x)\ Jb r (x) 

where fB R (x) s t an ds for the average of / over Br(x). 
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